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QPO as the Rosetta Stone for understanding black hole accretion 
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Abstract. Quasi periodic oscillations (QPO) seen in the X-ray fluxes of individual neutron stars and black hole sources are 
one of most intriguing phenomena in today's astrophysics. The QPO nature is visibly determined by super-strong Einstein's 
gravity. I argue here that it also profoundly depends on the MRI turbulence in accretion flows. Understanding the QPO 
physics may therefore guide accretion theory out of its present state of confusion. 
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1. Accretion theory in search of a paradigm 

Accretion provides energy for radiation in many astrophysi- 
cal sources: stellar binaries, active galactic nuclei, and proto- 
planetary disks. The black hole accretion in quasars is most 
powerful and most efficient steady energy source known in 
the Universe. In low-mass X-ray binaries (LMXRBs) the 
compact object, i.e. a neutron star (NS) or a black hole (BH) 
accretes mass lost by its binary companion, i.e. a normal, low 
mass star, M ~ M Q . BH sources in LMXRB are in many 
respects a scaled-down version of quasars, and they are often 
called microquasars 1 . 

Accretion theory experiences today a period of confu- 
sion: a picture of violently unsteady accretion on all its length 
and time scales emerging from supercomputer simulations, 
is plainly uncompatible with the long-lasting, quasi-steady 
accretion states (Shakura-Sunyaev, slim, adaf) predicted by 
semi-analytic methods. 

1.1. Shakura-Sunyaev disks, slim disks and adafs 

According to the semi-analytic understanding, developed 
over the past thirty years, the initially high angular momen- 
tum (of matter lost by the low mass star in an LMXRB) is 
gradually removed by viscous stresses that operate against 
the shear in the flow. This allows matter to gradually spiral 

Correspondence to: marek@fy.chalmers.se 
1 See an excellent non-technical review by Lasota(1999) 



down into the compact object, with gravitational energy de- 
graded to heat also by viscous stresses. A part of the heat 
converts into radiation, escapes, and cools down the flow. 

When radiative cooling is efficient, a Shakura-Sunyaev 2 
type of accretion flow is formed. It is optically thick but ge- 
ometrically thin (disk-like) and made of relatively cold gas 
(with negligible radiation pressure) that goes down on very 
tight spirals, resembling almost circular, almost free (Kep- 
lerian) orbits. Shakura-Sunyaev thin disks are relatively lu- 
minous, and have "thermal" electromagnetic spectra, i.e. not 
much different from that of a black body. 

Radiatively inefficient flows (RIF) are of two types — 
the first being a slim disk 3 , which is geometrically "slim" 
i.e. rather thin than thick, optically thick, and radiation pres- 
sure supported — and the second being an adaf 4 which is 
very hot, optically thin, and geometrically extended, similar 
in shape to sphere (or "coronae") rather than disk. Adaf emits 
a "power-law" non-thermal radiation, often with a strong 
Compton component. Slim disks are consistent with accre- 
tion rates of the order and higher than the Eddington rate, 
and Shakura-Sunyaev, and adaf are always sub-Eddington. 
Theory predicts that Shakura-Sunyaev, slim disks and adafs 
are stable: their characteristic steady states are much longer 
than the dynamical (Keplerian) times and comparable with 
thermal or viscous times at which a limit-cycle behavior may 

2 The classic work by Shakura & Sunyaev (1974) is one of the 
most often quoted papers in modern astrophysics. 

3 Abramowicz, Czerny, Lasota & Szuszkiewicz (1988) 

4 Abramowicz & al.; (1995); Narayan & Yi (1995). 
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1.2 The MRI instability 
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Fig. 1. The above "Bursa plot" illustrates a rather encouraging qualitative agreement between the twin peak QPO observations 
and the theoretical model based on a non-linear, strong-gravity resonance. Characteristically for a non-linear resonance the 
observed ratios of frequencies change in time but stay near the 3/2 value, determined directly by the strong gravity. For NS 
sources the frequencies locate approximately along the line v up = ao v\ + bo. This behavior is generic, and predicted by 
exact mathematics of the model. The strong-gravity 1/M scaling is also clearly apparent — see three black horizontal error 
bars that show measured masses of three microquasars (related to the upper horizontal mass axis). The NS QPO frequencies 
are about ten times higher than the BH frequencies, because NS masses are smaller by about the same factor. 



occur. For the same (sub-Eddington) accretion rate, Shakura- 
Sunyaev and adaf may simultaneously coexist at different 
radii or even overlap at some, and the flow may transient be- 
tween these two types due to some external instabilities (e.g. 
evaporation). In addition to the accretion inflow, outflows are 
present, in microquasars often in the form of transient, well 
collimated jets that shot out with speeds approaching that 
of light. Jets acceleration and collimation are, most proba- 
bly, governed by magnetic fields, but relevant details are un- 
known. It may be that the microquasars jets are powered by 
the Blandford-Znajek (1977) mechanism, an electromagnetic 
version of the Penrose process that utilizes BH's rotational 
energy (see also Li & Paczyriski, 2003). 



1.2. The MRI instability 

In a weakly magnetized fluid, with magnetic field B z and the 
Alvfen speed v\ = B 2 / (47rp) the dispersion relation for per- 
turbations of the fluid quantities SX ~ exp[i(fcz — cot)] reads 
cj 4 — (2kvA+^f) Lj2 + f tl 'A(kvA + rdfl 2 /dr) = 0, where fl is 
the angular velocity in the fluid, and u r is the radial epicyclic 
frequency. Obviously, this equation has an unstable solution 
u 2 < 0, if and only if, kvA + rdQ 2 /dr < 0. This is the 
celebrated magneto rotational instability (MRI), introduced 
to physics of accretion flows by Balbus & Hawley (1991; 
review 1998). MRI generates turbulence. Turbulent stresses 
dissipate orbital energy and redistribute angular momentum, 
thus allowing accretion to occur — without MRI there would 
be no accretion! However, in my opinion, this fundamental 
progress in understanding "from first principles" the nature 
of stresses in accretion flows is unmatched by present-day 
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supercomputer magnetohydrodynamical (MHD) simulations. 
In these simulations, accretion flows appear to violently vary 
at all time and length scales, never approaching a steady state. 

1.3. The need for a new paradigm 

The picture of violent accretion is obviously in a direct con- 
flict with the Shakura-Sunyaev, slim disk and adaf mod- 
els that predict long lasting and orderly quasi-steady states. 
There is no consensus on what does this acute rift between 
semi-analytic models and numerical simulations imply for 
our understanding of accretion. 

I agree with these who are convinced that the picture of 
quasi-steady accretion must be generally correct, and who 
point out several shortcomings of numerical simulations to- 
day. Simulations still depend too much on technical details 
and methods used (types of numerical boundary conditions, 
grid size), and results obtained by different groups do not 
agree in several relevant details. The simulations do not in- 
clude radiative processes, and make no predictions about 
spectral states of accretion flows. Despite their potential and 
trully impressive sophistication, they have not yet matured, 
and a lot of work is needed before they will. 

These who believe in MHD simulations, point out that 
while semi-analytic models will always depend on ad hoc 
phenomenological parameters to describe stresses, the MHD 
simulations calculate stresses directly from well understood 
piece of fundamental physics of MRJ. 

All in all, understanding of accretion is today somehow 
confused. Theory makes no sufficiently precise predictions 
to allow sharp observational judgements, and observations 
do not combine into a clear qualitative scheme. Phenomenol- 
ogy of spectral states that are really observed in LMXRBs 
is painfully complex. Neither it agrees with the "analytic 
paradigm" based on the three standard components (Shakura- 
Sunyaev — adaf or corona — jet), nor there is any indication 
that it could be explained in foreseeable future 5 by the "nu- 
merical MRJ paradigm". 

Major puzzles here are: what are the spectral states, what 
causes transitions between states, how transitions across the 
"jet line" help to launch a jet, and what is the physical reason 
for the connection between QPO and spectral states. 

In my opinion the last question is the key one. I am con- 
vinced that an answer to it would became the Rosetta Stone 
for accretion theory. 

2. The 3/2 twin peak QPO in LMXRBs 

Several of the observational facts discussed in this Section 
are illustrated in Figure ^ designed by Bursa (2003, unpub- 
lished). The Fourier power density spectra (PDS) of the ob- 
served LMXRBs X-ray variability reveal a broad band 1// 
noise, probably consistent with the MRI turbulence, and a 
few clear QPO peaks with quality factor Q > 10. Advanced 
data analysis points even to Q > 100. The highest observed 

5 J.-P. Lasota's estimate: "Maybe in ten years" 




Fig. 2. The epicyclic frequencies in Schwarzschild spacetime 

QPO frequencies are in the kHz range, corresponding to Kep- 
lerian orbital frequencies a few gravitational radii away from 
the central NS or BH. In many cases, "twin peak QPO" are 
present with two characteristic frequencies, v up and v\ Q , that 
have frequency ratio close to the v uv /v\ Q = 3/2 value 6 . The 
BH twin peak QPO have frequencies fairly fixed, with ratios 
being sharply at 3/2. The frequencies scale inversely with 
the BH mass 7 . The NS twin peak QPO frequencies wander 
in time in a wide range, but their changing ratios cluster re- 
markably sharply at 3/2, the same value as in the BH case. 
The NS data shows also the same 3/2 ratio in the anticorrela- 
tion (predicted by the resonance model) between slopes and 
shifts of the Bursa lines for individual NS sources. 

3. Epicyclic resonance model 

3.1. Observational motivations and simple physics 

The twin peak BH QPO are very accurate, high frequency, 
two coupled transient clocks, with frequencies fairly fixed 
for a particular source from one QPO event to another, that 
in addition scale from source to source as 1/M. In my opin- 
ion, any serious theoretician is forced to conclude from these 
facts that the QPO clock is determined by Keplerian motion 
in strong gravity, and not e.g. by MRI or some radiative pro- 
cesses. The simplest and the most natural possibility is that 
the clock mechanism is governed by small epicyclic devi- 
ations Sr, Sz from the exact circular (r — tq) and planar 
(z = 0) motion, as indeed expected in the Shakura-Sunyaev 
disk. In the zeroth-order approximation, the model thus ex- 
plains QPO in terms of two simple harmonic oscillators, con- 
sistent with the radial and vertical epicyclic motions, 

Sr + u% Sr = 0, Sz + u) 2 z 5z = Q, (1) 

where oj r and uj z are the radial and vertical epicyclic frequen- 
cies. In Newton's gravity with the spherically symmetric po- 

6 The 3/2 ratio for BH was first noticed, and its importance 
first stressed, by Abramowicz & Kluzniak (2001), and for NS 
by Abramowicz & al. (2003). Rational ratios were anticipated by 
Kluzniak & Abramowicz (2000). 

7 The 1/M scaling was first found by McClintock & Remillard 
(2003) for the three microquasars with known masses (see Fig.0- 
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tential it is u> r = uj z — lok, where lok is the Kepler's orbital 
frequency, but in strong gravity (see Figure^, 

U z > U K > UJ r , (2) 

and w r (rMs) =0, u 2 < for r < tms, where t"ms = ISCO is 
the radius of the marginally stable circular orbit. 

For a nearly-Keplerian^M/a! (Shakura-Sunyaev disk, slen- 
der torus), with the fluid sound speed C| /uj 2 r 2 = (3 -C 1, the 
frequencies of fluid's epicyclic oscillations Cu r , uj z are modi- 
fied by the pressure 8 , 

Q r = lu t - XrP, ui z =u z - XzP, Xr,Xz = const (3) 

The modification is very small, but crucial, as it provides 
a weak pressure coupling between the epicyclic modes of 
fluid oscillations that leads to the 3/2 resonance. To see 
how it works, one writes Q and |T5} in Mathieu's form, 
5'z + lu 2 [1 + xcos(Cu r t)]5z = 0, Sr + Co 2 5r = 0. A para- 
metric resonance occurs when uj r ~2uj z /n, n— 1, 2, 3, ... In 
strong gravity u> r <u) z , and thus n=3 is the lowest possible 
value of n, corresponding to the strongest resonance. This ex- 
plains why the ratio, uj z /uj r = ^up/^io = 3/2 is most often 
observed. 

3.2. Mathematical formulation of the model 

The Kluzniak & Abramowicz QPO resonance model starts 
from Vj T\ = 0. Physics, however complicated, is fully de- 
scribed by stress-energy tensor T l k (t) = °T\ + N l k {t). The 
"mean-field", statistical average of the flow is described in 
terms of a slowly varying tensor °T l k . Turbulent stresses (e.g. 
the MRI turbulence), fluctuate in time and are described by 
Af\(t). A similar division of T\ was employed, for a dif- 
ferent purpose, by Ogilvie (2003). Deviations from a quasi- 
Keplerian flow are described by <$V, T \ = 0, which is put 
into, 

5r + lj 2 Sr = F(Sr, Sz, Sr, Sz) + A cos(uj Q t) + J\f(t), 

Sz + lo 2 z Sz = G{Sr,Sz,Sr,Sz) + Bcos(io t) +JV(t). (4) 

Here the functions F, G account for coupling (known is terms 
of expansion of SVi°T\ = 0), the cos(o;ot) terms describe 
external forcing (relevant in the NS case, loq is the NS spin), 
and Af(t) describes an influence of the MRI turbulence. Ob- 
viously, cos(cjo^) an d N(t) f ee d energy into the resonance. 

Equation @ was studied by the method of multiple 
scales. All variables, including time, were expanded accord- 
ing to X = ^2 £ k Xk, e -c 1, and equations of different orders 
of e solved separately. Some of the coefficient in expansions 
have been assumed ad hoc, as it is difficult to derive all of 
them from first principles. 

The 3/2 resonance was found to be a natural and robust 
property of nearly Keplerian fluid motion is strong gravity, 
and resulting from it double peak QPO frequencies have been 
found to behave as these observed, in particular the Bursa line 
was recovered, see Figure^ 

Obviously, one wants derive all expansion coefficient 
from first principles. Of a particular interest is the order 

8 As first explained in the QPO context by Kluzniak & Abramow- 
icz (2002). Discussion that follows is copied from their paper. 



Af(t) ~ that should follow from a quantitative under- 
standing the MRI turbulence. A pure mathematical reasoning 
hints how a physically interesting range of N may be found. 
When N is too large the turbulence is too small and supplies 
not enough energy to keep the resonance alive. When N is 
too small, the turbulence is too strong and dominate dynam- 
ics, preventing QPO from showing up above the noise. Thus 
TV must be in a particular range, that may be determined e.g. 
by stochastic differential equations (SDE) that we use. One 
hopes to compare this range with that predicted by MRI tur- 
bulence. 

4. Forced oscillations: QPO and NS spin 

In in some NS sources, the QPO frequencies obviously de- 
pend on the NS spin. For example, difference in frequencies 
of the double peaked QPO in the millisecond pulsar SAX 
J1808. 4-3658 is clearly equal to half of the pulsar spin. This 
suggests that the epicyclic resonance is excited by direct forc- 
ing of accretion disk oscillation modes by the neutron star 
spin. 

We studied the forcing in the case of slender tori around 
rotating NS. All possible linear normal modes of oscillations 
of a perfect fluid slender torus are known analytically (Blaes, 
1985). They include the "epicyclic" fluid modes described by 
equation @, also in the case of a slightly non-slender torus. 
We study analytically non-linear behavior of these modes, in- 
cluding their coupling and resonances, MRI turbulence, and 
direct external forcing cos(ojot) by the NS rotation ujq. This 
is a rather difficult mathematical problem that takes time — 
first results are excepted this year 9 . So far, all our published 
results on the cos(wot) external forcing of oscillations of tori 
are based on numerical simulations. 

It was found that a resonant response occurs when the 
difference between frequencies of the two epicyclic modes 
equals to one-half of the spin frequency (as observed in SAX 
J1808. 4-3658 and other "fast rotators"), and when it equals to 
the spin frequency (as observed in "slow rotators" like XTE 
J 1807-294). 

5. Spectral states, QPO, modulation of X-rays 

Jean-Pierre Lasota told me recently, reflecting on the obser- 
ved puzzling connections between QPO and spectral states of 
LMXRB — "I am convinced that resonances contribute only 
a small part to something that is more fundamental there, and 
connected to how does the disk radiate... 10 ". Nobody knows 

9 Our unpublished results, both analytic and numerical, revealed 
a remarkably complex network of mode coupling and resonances 
that occur when forcing is included. A similar numerical research 
of oscillations of fluid tori is being conducted by Omer Blaes and 
Chris Fragile at UCSB. Rezzolla & al. (2003) were first to study 
oscillations of tori in the QPO context. 

10 Jestem przekonany, ze rezonanse sa tylko malym przyczynkiem 
do czegos bardziej podstawowego, zwiazanego z tym jak dysk 
swieci, a MRI (na razie) nie modeluja dyskow i nie mowia nic o 
ich swieceniu. Moze za dziesiec lat symulacje beda przedstawialy 
dyski. Takie jest tez zdanie Boba W. 
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how to explain the puzzle, and we are all guessing. My in- 
tuition tells that in both NS and BH sources the QPO clock 
must be regulated by accretion disk oscillations that are set by 
strong-gravity's orbital dynamics. In my opinion, how does 
the matter radiate cannot be fundamentally important for the 
QPO clock. However, the observational appearance of QPO 
must be influenced by the same processes that modulate radi- 
ation, and this is where I could see a natural link between 
QPO and spectral states. To be more specific, let me start 
from explaining how some differences in appearance of QPO 
in BH and NS sources may follow from the "unique gravity's 
clock — non-unique modulation" principle. Close to ISCO, 
the vertical oscillation 5z(t) modulate Mi n which is very sen- 
sitive to the mariginally stable ISCO conditions, 

(5M/M) m = X Sz 2 , x = X(r%n,n) = const. (5) 

In NS sources, the total X-ray luminosity Lx is dominated 
by that of the boundary layer, Lx ~ L X L ~ M$ n . Thus, the 
QPO modulation of the luminosity occurs at the boundary 
layer. 

For BH this mechanism does not work because the mod- 
ulated inner mass flux M- m disappears inside the BH event 
horizon, and thus no (or only a very little) modulation of the 
X-ray luminosity may be expected. However, it was shown 
by a large scale ray tracing method that in this case, the X- 
ray luminosity is modulated by strong gravity effects (gravi- 
tational lensing and Doppler beaming) in photon propagation 
from oscillating disk to observer. It was found that the "ob- 
served" power in PDS Fourier spectra very strongly depends 
on lensing and beaming effects — for example, power in the 
vertical mode increases with inclination. 

The ray tracing simulations proved that observed PDS 
power of QPO modes depends also on the vertical size of 
the oscillating region of the disk. Obviously, the size varies 
with temperature that is an important parameter for spectral 
states. This provides an (unavoidable, based only on photon 
propagation in strong gravity) link between spectral states 
and QPO. Note, that it is believed that in different spec- 
tral states the location of inner radius of the optically thick 
Shakura-Sunyaev disk varies. This provides an opaque vary- 
ing screen that blocks some of the photons involved in lens- 
ing and beaming, with important consequences for the QPO 
visibility in different spectral states. 

Of course, QPO appearance may also be intrinsically in- 
fluenced by different physical conditions of the flow in dif- 
ferent spectral states. This we study by a direct ray tracing 
from Mami Machida's accretion flow MHD 3D simulations 
in Schwarzschild spacetime. 

6. Conclusions 

Similarities in QPO observed in BH, NS and WD sources 
point to a common origin — accretion disk oscillations. Dif- 
ferences reflect different masses and sizes of BH, NS and WD 
that induce differences in excitation of oscillations and mod- 
ulation of luminosity. 

In particular, repeatability and precision of the two QPO 
clocks in double peak QPO in BH and NS sources is ex- 



plained by a resonance between two epicyclic modes, ver- 
tical and radial, that exist in all nearly Keplerian fluid accre- 
tion flows. The 3/2 frequency ratio and the 1/M scaling are 
both a direct and natural (unavoidable) consequence of the 
strong gravity. The resonance is fed by MRI turbulence in 
BH disks, and may be additionally forced by rotation of the 
central star in NS sources. Modulation of X-ray luminosity 
occurs at the boundary layer in NS sources, and is due to 
lensing and Doppler beaming in BH sources. Observed cor- 
relations between QPO and the LMXRB spectral states do 
not necessarily imply that QPO strongly depend on radiative 
processes — more likely the correlations are determined by 
the QPO "visibility" during different spectral states. 

Nature and properties of double peak QPO in BH and 
NS sources, profoundly depend on Einstein's strong gravity, 
and on local physics of the MRI turbulence. For this reason, 
a deeper understanding of QPO may provide a rescue from 
confusions that upset the accretion theory today. 
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